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Abstract 

In recent years, single molecule force techniques have opened a new avenue to decipher the folding 
landscapes of biopolymers by allowing us to watch and manipulate the dynamics of individual pro- 
teins and nucleic acids. In single molecule force experiments, quantitative analyses of measurements 
employing sound theoretical models and molecular simulations play central role more than any other 
field. With a brief description of basic theories for force mechanics and molecular simulation technique 
using self-organized polymer (SOP) model, this chapter will discuss various issues in single molecule 
force spectroscopy (SMFS) experiments, which include pulling speed dependent unfolding pathway, 
measurement of energy landscape roughness, the influence of molecular handles in optical tweezers on 
measurement and molecular motion, and folding dynamics of biopolymers under force quench condition. 
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The advent of single molecule (SM) techniques over the past decades has brought a significant 
impact on the studies of biological systems [IH3]. The spatial and temporal resolutions and a 
good force control attained in SM techniques have been used to decipher the microscopic basis 
of self-assembly processes in biology. Among SM techniques, single molecule force spectroscopy 
(SMFS) have been adapted not only to stretch biopolymers [HE] but also to unravel the internal 
structures and functions of many proteins and nucleic acids [dTTIU]. By precisely restricting the 
initial and final conformations onto specific regions of energy landscape, SMFS also provided a 
way to probe the collapse or folding dynamics under mechanical control, which fundamentally 
differs from those under temperature or denaturant control [TT] [12] . The observables that are 
usually inaccessible to conventional bulk experiments, for instance the heterogeneity of dynamic 
trajectories and intermediate state ensembles, have been measured to provide glimpses to the 
topography of complex folding landscapes [13]. Furthermore, the use of SMFS is being 
expanded to study the function of biological motors [HH2Q] and cells [2T] . 

Given that foldings of biopolymers are realized through a number of elementary processes, 
good controls over time, length, force and energy scales are essential to resolve the details of 
biomolecular self-assembly [22H21]- The ability to control the energy scale within the range 
of ~ ksT, in particular, (ksT ^4.1 pN ■ nm at room temperature T = 300 K) allows us 
to study how biological systems, that are evolved to accommodate the thermal fluctuations, 
versatilely adapt their conformation to a varying environment. SMFS is an excellent tool to 
decompose the energy required to disrupt non-covalent bonds, responsible for the stability of 
biological structures (~ 0(1) /c^T), into ~pN force and ~nm length scale. A phenomenological 
interpretation of bond rupture due to an external force in the context of cell-cell adhesion 
process and theoretical estimate of mechanical force associated with the process had already 
been discussed as early as in 1978 by Bell |25j. However, only after 1990s with SMFS its 
experimental realization was achieved [26-30J. 

Many biological processes, in vivo, are in fact mechanically controlled. The ability to apply 
pN-force to a single molecule and watch its motion at nm scale (or vice versa) has a great 
significance in molecular biology in that one can elucidate the microscopic and structural origin 
of a biological process by quantify both kinetic and thermodynamic properties of the biomolecule 
at single molecule level and compare them with those from ensemble measurements [31] . Under 
a constant force condition using force-clamp method [11] HH [231 E2] near at the transition region, 
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the molecular extension (x) exhibits discrete jumps among basins of attractions as a function of 
time. This immediately allows to study the hopping kinetics between the basins of attraction. 
If the time traces are long enough to sample all the conformations then one can also construct 
an equilibrium free energy profile under tension / [331 US]- The fraction of nativeness 0jv(/) or 
equilibrium constant K eq (f) = (1 — 0at(/))/0at(/) as a function of / can be accurately measured 
just like the one using calorimetry or denaturant titration in bulk experiments. 

Since the first single molecule force experiment, interplay between theory, simulation and 
experiment have affected the experimental design as well as theoretical formulation to interpret 
results from measurements. A need to understand biomolecular dynamics at SM level further 
highlights the importance of theoretical background such as polymer physics [35J , stochastic the- 
ory [SHI EZ] and fluid dynamics. Molecular simulations of SMFS using a simple model provide 
a number of microscopic insights that cannot be easily gained through experiments alone. This 
chapter encompasses the force mechanics from the perspectives of theories and molecular simula- 
tions. Basic theories for force mechanics and the main simulation technique using self-organized 
polymer (SOP) model will be described, followed by a number of findings and predictions for 
SMFS made through the concerted efforts using force theories and molecular simulations. 

Force Extension Curve 

Mechanical response of a molecule is expressed with two conjugate variables, force (/) and 
molecular extension (x), to define a mechanical work (W = f ■ x). Force-extension curves 
(FECs) or the time dependences of f(t) and x(t) are the lowest level data that can infer all 
the relevant information concerning the mechanical response of the biomolecules. For a generic 
homopolymer whose Flory radius is Rp ~ N u a, where N is the number of monomers and a is 
the size of monomer, the extension of polymer x should be determined by a comparison between 
Rf and tensile screening length £ p = ksT/f. The force value / determines the parameter 
q = Rp/^p, which satisfies q <C 1 (q 1) for small (large) /. The applied force can roughly be 
classified into three regimes, (i) For small /, x <C Rf and q are satisfied. Thus, x ~ (3R F f 
is obtained from a scaling argument x = Rp<&{q) ~ RfQ since 3>(q) ~ q for q -C 1. (ii) For an 
intermediate / (Rf < x <C Na), the shape of globular polymer is distorted to form a string of 
tensile blobs, where blob size is given £ p ~ N£ with being the number of monomers consisting 
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the blob. The total extension of the string of tensile blobs under tension is x ~ £ p x (N/N},), 
leading to Pincus scaling law x ~ Z 1 ^ -1 ~ / 2 ^ 3 [38J. Here, note that this scaling law is only 
observed when 1 <C (Cp/ a Y^ u ^ N is ensured [39j. (iii) For extremely large forces, chain is fully 
stretched; x ~ Na 2 (3f /3 for an extensible chain and x ~ Na for an inextensible chain. 

The scaling argument for biopolymers deviates from that of generic homopolymers with 
N ^> 1 due to the finite size effect as well as various local and nonlocal interactions [39]. In 
practice, the persistence length (l p ) and contour length (L) of biopolymers are extracted by 
employing a force-extension relation (/ = [1/4(1 — x/L) 2 — 1/4 + x/L] [10]) of worm-like 
chain (WLC) model, whose energy hamiltonian takes into account the bending energy penalty 
along the polymer chain (H/kgT = ^ J Q L ^ 9 ^ j ds where u(s) is the tangential unit vector 
at position s along the contour) [I]. The rips in FEC due to the disruption of internal bonds 
and subsequent increase in the contour length from L to L + AL are used to decipher the 
energetics and internal structure of proteins and nucleic acids. [6j [HI [23] . The FEC of repeat 
proteins demonstrates multiple peaks with saw-tooth pattern, suggesting that under tension the 
repeat proteins unfolds one domain after another. As a more complicated system, Tetrahymena 
ribozyme with nine subdomains show saw-tooth patterns but with varying peak height and 
position, demanding more careful and laborious tasks of analysis [5]. 

Forced Unfolding at Constant Force 

A phenomenological description of the forced-unbinding of adhesive contacts by Bell [22] has 
played a central role in studying the force induced dynamics of biomolecules for the last two 
decades. In the presence of external force /, Bell modified Eyring's transition state theory jH] 
as follows : 

where ks is the Boltzmann constant, T is the temperature, h is the Planck constant, and k is 
the transmission coefficient. The parameter 7 is a characteristic length of the system associated 
with bond disruption. Under tension /, the activation barrier responsible for a stable bond is 
reduced to Eft — 7 X /. The prefactor is the vibrational frequency of a bond due to thermal 
fluctuation prior to disruption. 

Although the original Bell model correctly describes the stochastic nature of bond disrup- 
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tion, the prefactor ksT/h fails to capture the physical nature of attempt frequency for the 
ligand unbinding from catalytic site or the protein unfolding dynamics under tension, which 
depends on the shape of potential as well as viscosity of media. In fact, Eyring's transition 
state theory is only applicable for the chemical reaction in gas phase. More appropriate the- 
ory for the dynamics in condensed media should account for the effect of solvent viscosity and 
conformational diffusion [32]. One-dimensional reaction coordinate, projected from a multidi- 
mensional energy landscape, can well represent the dynamics provided that the relaxation times 
of conformational dynamics along a reaction coordinate is much slower than other degrees of 
freedom [13]. Under tension /, molecular extension (x) (or end-to-end distance (R)) are as- 
sumed to be a good reaction coordinate. On the one-dimensional reaction coordinate, mean 
first passage time obeys the following simple differential equation [33]. C fp (x)t(x) = — 1 where 
C FP = e Feff( - x ^ kBT d x D(x)e~ Feff( - x ^ hBT d x is the adjoint Fokker Planck operator. The mean first 
passage time of a quasi-particle between the interval a < x < b with reflecting d x r(a) = and 
absorbing boundary condition r{b) = is 



Above, the free energy profile F(x) is considered being "tilted" by an external force by / -x. As 
long as the transition barrier (AF* = F(x ts ) — F(xf,)) is large enough, the Taylor expansions 
of the free energy potential F(x) — fx at the barrier top and the bound state position with a 
saddle point approximation result in the seminal Bell-Kramers equation [3SI EZ] , 



where Aa;* = x ts — Xb, ujb and u ts are the curvatures of the potential, \d%F(x)\, at x = x^ 
and Xts, respectively, and 7 = ksT/Dm is a friction coefficient associated with the motion of 
biomolecule. Experimentally determined speed limit of the folding dynamics (barrierless folding 
time) of two-state proteins is ~ (0.1 — 1) fis [321 133 EE]. A care should be taken not to use 
the prefactor (kpT/h)^ 1 0.2 ps from transition state theory for gas phase when estimating 
the barrier height from folding or unfolding kinetics data of biopolymers in condensed phase. A 
cautionary word is in place. If the barrier height AF* — /Aor is comparable or smaller than 
ksT, the molecular configuration trapped as a metastable state in the free energy barrier can 
move almost freely across the barrier. In face, the barrier vanishes when / reaches a critical force 
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f e = AF*/Ax*. In this case, there is no separation of time scales between diffusive motion and 
barrier crossing event. Hence, the saddle point approximation taken for Eq|3]from Eqj2] does 
not hold. Due to thermal noise the unfolding or rupture event of the system occurs under finite 
free energy barrier (> fc^T) before / reaches f c . For Bell-Kramers equation to be applicable / 
should be always smaller than f c . 



Forced Unfolding at Constant Loading Rate - Dynamic Force Spectroscopy 

Even though the constant force (force clamp) experiment is more straightforward for analysis, 
due to technical reasons many of the force experiments have been performed under a constant 
loading condition (force-ramp) in which the force is linearly ramped over time P, H7H3H]. Dy- 
namic force spectroscopy (DFS) probes the energy landscape of biomolecular complexes by de- 
tecting the mechanical response of the molecules. The linearly increasing mechanical force with 
a rate of rf = df/dt is exerted on the molecular system until the molecular complex disrupts. 
Upon unbinding, the force recorded on the instrument drops abruptly, thus one can measure the 
unbinding force of the system of interest (Fig,lA). Because of stochastic nature of unbinding 
event, the unbinding force of molecular complex is not unique, but distributes broadly, defining 
the unbinding force distribution (P(f)) (Fig. IB). 

Under linearly varying force (/ = rj x t), the rate of barrier crossing from bound to unbound 
state (or from folded to unfolded state) is also time-dependent. Hence, the survival probability 
at time t is given by S(t) = exp (— j^drk(T)\. Thus, the first passage time distribution is 
P(t) = —dS(t)/dt = k(t)S(t). Change of variable from t to / leads to a unbinding force 
distribution 

i i r rf i i 

(4) 



P(f) = —k(f)S(f) = —k(f) exp 
r f r f 



f 1 

/ df-k(f') 
Jo r f 



Note that k(f) is exponentially increasing function of / while S(f) is exponentially decreasing 
function of / with greater power at / 3> 1, shaping a Gumbel distribution, P(f) ~ e^e~ ef , 
for k(f) ~ e* . Current theoretical issue of deciphering the underlying energy landscape using 
force hinges on an analysis of P(f) by building not only a physically reasonable but also a 
mathematically tractable model. 
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f = A^ l0gr/+ A^ l0g ^ 



(5) 



The most probable unfolding force is obtained using dP(f)/df\f=f* = 0. 

Ax t 
u D e-^ Ft k B T 

where = ^o^ts/^l- in conventional DFS theory, Eq|5] is employed to extract Ax$ and 
AF* of underlying 1-D free energy profile associated with force dynamics. For unbinding force 
/* to be compatible to the one in the picture of Kramers barrier crossing dynamics, /* < f c 
should be obeyed as mentioned above. The condition /* — f c — log T J^ T < demands 
Tf < rj (= UD ^ t T ) ■ For a set of parameters, k&T w 4 pN ■ nm, Ax^ ~ 1 nm, and vq ~ 10 6 
the critical loading rate is ~ 10 6 pN/s. The typical loading rate used in force experiments 
(O.lpN/s < Tf < 10 3 pN/s) is several orders of magnitude smaller than this value. Therefore, 
in all likelihood unbinding dynamics in typical experimental conditions obey stochastic barrier 
crossing dynamics. In contrast, the steered molecular dynamics simulations with all atom rep- 
resentation [50J typically uses r/ > due to high computational cost. In such an extreme 
condition, however, the forced-unfolding process can no longer be considered a thermally acti- 
vated barrier crossing process. At high ry > it was shown that an average rupture force ((/)) 
grows as [5T] . 

It is of particular interest that for a molecular system unfolding through a single free energy 
barrier, the force dependence of force clamp kinetics can be formally expressed with the P(f) 
from force ramp experiment as follows 



k(f)= min (Q) 

l-fo f df>P(fi)/r f 

which is easily shown using the relation, S(f) = exp (— Jq df'^jk(f') \ = 1 — Jq df'^P(f'). 
Technically the two distinct experimental methods are connected through this simple relation- 
ship. Therefore, by conducting force-ramp experiment with a sufficiently good statistics to get 
P(f) at varying r/, one can, in principle, build a data for k(f) as in force-clamp experiment. 



Deformation of Energy Landscape under Tension 



Basic assumption of Bell-Kramers equation is that an external force changes the free energy 
barrier along the reaction coordinate from AF$ to AF$ — /Aar without significantly chang- 
ing other topology of energy landscape. A linear regression of both Eqs. [3] and [5] provides 
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the characteristic length Ax* and free energy barrier AF* at / — > 0. However, in practice, 
nonlinearity (negative curvature for / vs log k(f), positive curvature for log 77 vs /*) is often 
detected especially when / or r/ is varied over broad range. Thus, in case / (or r/) is varied 
at a large but narrow range of / (or rf) then substantial errors can arise in the extrapolated 
values of Ax* and AF* to the zero force ; the linear regression will underestimate Ax* and 
overestimate AF*. The physical origin of /-dependent Ax* is found in a complicated molec- 
ular response to the external force. If the TS ensemble is broadly spread along the reaction 
coordinate then the molecule can adopt diverse structures along the energy barrier with varying 
/ values. Whereas, if the TS ensemble is sharply localized along the reaction coordinate, the 
nature of TS ensemble measured in x-coordinate will be insensitive to the varying / values. Or, 
more simply, the origin of moving transition state position can be algebraically explained by 
plotting the shape of F e ff(x) with varying /. Because Xt s and x b are determined from the force 
dependent condition F'(x) — f = 0, all the parameters should be intrinsically /-dependent as 
AF*(/), Ax*(/), Uts{f), and 0Jb{f)- By making harmonic approximation of F(x) at x = Xb and 



x = x tB , i.e., F(x) « F(x b ) + 1/2 • \F"{x b ) | (x - x b f and F(x) « F(x ts ) - 1/2 • \F"(x ts ) | (x - x ts f 



and calculating Ax*(/) = x ts (f) — x b (f) from F^Jx) = F'(x) — / = 0, one can show that 



Typically for biopolymers under tension, free energy profile near native state minimum is sharp 
(\F"(x b ) \ 3> 1). To minimize the difference between Ax*(/) and Ax* and to make x(f) ~ 0> the 
transition barrier should be sharp for a given //Ax* value (i.e. //Ax*|F"(x£ S )| <C 1). Simulation 
studies [55H57] in which the free energy profiles were explicitly computed from thermodynamic 
considerations alone clearly showed the change of Ax* when / is varied. Note that the movement 
of transition barrier location toward the native state position (x = x b ) is consistent with the 
Hammond postulate EH] that explicates the nature of transition state of a simple organic 
compound when product state is relatively more stabilized than reactant state. 

To account for the nonlinear response of biological systems to the force more naturally, Dudko 
and coworkers proposed to use an analytically tractable microscopic model for the underlying 
free energy profile. For a cubic potential F(x) = — x 2 (x — |Ax*) whose distance to the 
transition state and free energy barrier are Ax* and AF*, all the parameters needed for Kramers 
equation are expressed as a function of /, Ax* and AF*; u> ts (f)uj b (f) = u ts Ub(l — f/fc) 1 ^ 2 , 
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Ax*(f) = Az*(l -f/f c ) 1/2 , and AF i (f) = |A^/ C (1 - /// c ) 3/2 where / c = 3AF*/2Aa;t. Thus, 
the /-dependent unfolding rate k(f) for cubic potential is given exactly with v = 2/3. 



k(f) = k(0) ( 1 - J-^ ) e 



AFT 



In fact, with different v value, the same expression with Eq|8] is obtained for harmonic-cusp 
potential [y = 1/2) and Bell-Kramers model [y = 1). The precise value of v depends on 
the nature of the underlying potential and could be treated as an adjustable parameter [53J. 
Consequently, cusp, cubic [52j |53j [60] or piecewise harmonic potentials [61] have been suggested 
as microscopic models for underlying free energy profile. So far, theories for force experiments 
have been devised mainly for single barrier picture in one-dimensional reaction coordinate. 
Although two slope fit using multiple energy barrier picture was suggested to explain the large 
curvature observed in DFS data [62J, building an multibarrier free energy profile from one- 
dimensional information such as P(f) or [logr/,/*] curve is an inverse problem whose answer 
may not be unique, as is well demonstrated by Derenyi et al. in the context of the forced 
unfolding over two sequentially located transition barriers 



Self-Organized Polymer (SOP) Model for Single Molecule Force Spectroscopy 

At present, the typical spatial resolution reached in SMFS is a few nm, and the dynamics 
that are probed in SMFS are rather global than local. The details of local dynamics such as the 
breakage of a particular hydrogen bond or an isomerization of dihedral angle cannot be discerned 
from FECs or the time trace of molecular extension alone. To gain and provide sufficient insight 
through molecular simulations in conjunction with SMFS, a simple model with which one can 
efficiently simulate the forced unfolding dynamics of a large biopolymer at a spatial resolution 
of SMFS would be of great use. By drastically simplifying the details of local interactions such 
as bond angle or dihedral angle along the backbone, which indeed are the major determinant for 
the dynamics under tension at nm scale, one can either gain an acceleration in simulation speed 
or explore the dynamics of a larger molecule. As is well appreciated in the literature of normal- 
mode analysis, an inclusion of small length scale information do not alter the global dynamics 
corresponding to a low frequency mode [5H45B] . Self-organized polymer model, proposed in this 
line of thought, is well suited to simulate the forced unfolding dynamics of large biopolymers 
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at spatial resolution of SMFS. The basic idea of the SOP model is to use the simplest possible 
energy hamiltonian to faithfully reproduce the topology of native fold and to simulate the low- 
resolution global dynamics of biopolymers of arbitrary size [67H73] . The energy function for 
biopolymers in the SOP representation is 

EsoPm) = E FENE + E£» + 

2^o log(l M ) 

i=i u 



N-3 N 
i=l j=i+3 



Kj 

The first term in Eq. [9] is the finite extensible nonlinear elastic (FENE) potential for chain 
connectivity with parameters, k = 20kcal / (mol ■ A2), Rq = 0.2 nm, r^+i is the distance between 
neighboring beads at i and i + and rf i+1 is the distance in the native structure. The use of the 
FENE potential for backbone connectivity is more advantageous than the standard harmonic 
potential, especially for forced- stretching to produce an inextensible behavior of WLC The 
Lennard- Jones potential is used to account for interactions that stabilize the native topology. 
A native contact is defined for bead pairs i and j such that \i — j\ > 2 and whose distance is 
less than 8 A in the native state. We use = 1 — 2 kcal/mol for native pairs, and e% — 1 
kcal/mol for nonnative pairs. In the current version, we have neglected nonnative attractions. 
This should not qualitatively affect the results, because under tension such interactions are 
greatly destabilized. To ensure noncrossing of the chain, + 2) pairs interacted repulsively 
with a = 3.8 A. There are five parameters in the SOP force field. In principle, the ratio of 
€h/ei and R c can be adjusted to obtain realistic values of critical forces. For simplicity, we 
choose a uniform value of th for all protein constructs, can be made sequence-dependent and 
ion-implicit as — > if one wants to improve the simulation results. By truncating forces 
due to the Lennard- Jones potential for interacting pairs with r > 3r° or 3cr, the computational 
cost essentially scales as ~ O(N). We refer to the model as the self-organized polymer (SOP) 
model because it only uses the polymeric nature of the biomolecules and the crucial topological 
constraints that arise from the specific fold. 

With SOP representation, the dynamics of biopolymers are simulated under force-clamp or 
force-ramp condition by solving the equation of motions in the overdamped regime. The position 
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of the i th bead at time t + h is given by 



Xi {t + h)= Xi (t) + - (Fi(t) + Ti(t) + f N (t)) 



(10) 



where Fi(t) = — V Xi £^5 , op({^}) ) the x-component of the conformational force acting on the 
bead, and T(t) is a random force selected from a Gaussian noise distribution P[Tj(t)] oc 



;th 



exp 



i 



4k B T( 



T 



. The conversion of simulation time into the physical time is made by 
using d t r sia/r~ C _1 (/c B T/a) and ma 2 /eh = r|, 

(a 2 \C{r L /m)-e h 



k B T 



k B T 



Ti- 



ll) 



When amino acid residue is used as a coarse-grained center, ( (50 — 100) tu/tl. The time step 
h ~ 0.1r£. For force-clamp condition, a constant force (/at) is exerted to the N th bead fisr(t) = 
/jv with the 1st bead being fixed, f\{t) = 0. For force-ramp condition, /jv(t) = —k(xN — vt) 
with fi(t) = is used. A harmonic spring with stiffness k is attached to the N th bead and the 
position of spring is moved with a constant velocity v. 



Deciphering the Energy Landscape of Complex Biomolecules using SMFS 

The mechanical response of the molecule becomes more complex with an increasing complex- 
ity in the native topology of biomolecules. This section will address two classes of unfolding 
scenarios for the molcules with complex topology. 

First, it is conceivable that a biomolecule with complicated topology in its native state is 
unravelled via more than two distinct transition state ensemble so that the forced-unfolding 
routes bifurcate [7J [71] . In the scenario, the survival probability of molecule remaining in the 
native state decays as 

N 

S N (t) = J2^(f)exp(-t/r l (f)) (12) 

i=l 

where unfolding time along the z-th route is given by Tj(/) and tpi(f) is the partition factor for 
z-th route with YliLi = 1> both of which are the function of / [12"1 [75] . Mechanical response 
of a barrel shaped Green fluorescence protein (GFP), made of 11 /3-strands with one a- helix at 
the N-terminal, is quite intricate, whose unfolding path depends on pullling speed and direction 
[7J [76]. In earlier force experiment on GFP by Rief and coworkers, the intricacy of the GFP 
forced-unfolding is manifested as the indistinguishability of two unfolding routes [77] . After the 
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a-helix being disrupted from the barrel structure due to external force, the second rip is due 
to the peeling off of (31 or /311. The gains of contour length (AL) from the molecular rupture 
event from (31 and from 011 are, however, identical, so that it was impossible to identify the 
source of the second peak simply by analyzing the FEC [77]. The force simulation using SOP 
model suggested a bifurcation into the two different routes by showing that 70 % of the molecule 
disrupt from N-terminal and the remainder of the molecules from C-terminal. Experimentally, 
this is confirmed by fixing either N-terminal or C-terminal direction by introducing a disulfide 
bridge through mutation [7J. The kinetic partitioning mechanism (KPM) used for protein and 
RNA folding [121 [78] can be adapted to explain the mechanical behavior of biopolymers. 

Second, reaction coordinate under tension can have sequentially aligned multiple barriers. 
In comparison to GFP, the forced-unfolding experiment on RNase-H was characterized by a 
peculiar mechanical response [79]. The FEC of RNase-H has a single large rip along unfolding 
path but has two rips in refolding path. This behavior was explained by considering a shape of 
free energy profile where native (N), intermediate (J), and unfolding (U) state lie sequentially 
with relatively high transition barrier ksT) between N and I, and low transition barrier 
(~ IvbT) between / and U. On such a free energy profile, the unfolding to the U state would 
occur by skipping I state because the external force that disrupt to the iV to overcome the 
first barrier is already larger than the mechanical stability of the I state relative to U. Thus, 
Accessing to the I state is difficult from the N under an increasing tension while I can be 
reached from U in the refolding FEC since the free energy barrier between I and U is relatively 
small (~ ksT). This hypothetical picture was further supported by the force-clamp method. 
For I ^ U, the transition mid-force is f m ,i^u ~ 5.5 pN while the escape from native basin of 
attraction (NBA) occurs at / « 15 — 20 pN |79j. Even for a molecule hopping between / and 
U at / = 5.5 pN, the molecule can get to N state. However, once the molecule jumps over the 
barrier between N and I, being trapped in N. the molecule has little chance to jump back to I 
within the measurement time. 

Measurement of Energy Landscape Roughness 

Although the energy landscape of biopolymers are evolutionary tailored such that the poten- 
tial gradient toward the native state is large enough to drive the biopolymers to their native state, 
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the energetic and topological frustration still remain to render the folding landscape rugged, slow- 
ing down the folding processes. To account for the effect of energy landscape roughness in a one- 
dimensional free energy profile, F(x) can be effectively decomposed into F(x) = Fq(x) + Fx(x) 
[80] . where F (x) is a smooth potential that determines the global shape of the energy landscape, 
and Fi(x) is the ruggedness that superimposes F (x). By taking the spatial average over Fi(x) 
using (e ±l3Fl ^}i = j J dxe ± P Fl ( x \ where I is the ruggedness length scale, the associated mean 
first passage time is altered to r(x) w j b x dye Fo{y)/kBT {e m{y) )i^ dze- Fo ^/ kl}T {e~ pFl ^)i. By 
either assuming a Gaussian distribution of the roughness contribution F\ (P(Fi) oc e~ F ^ 2e2 ) 
or simply assuming (3F\ <C 1 and (Fi) = 0, (F 2 ) = e 2 and (3e is small, the effective diffusion 
coefficient can be approximated as D* ~ Dexp (— (3 2 e 2 ) where D is the bare diffusion constant. 

The signature of the roughness of the underlying energy landscape is uniquely reflected in 
the non-Arrhenius temperature dependence of the unbinding rates. By using mean first passage 
time with the effective diffusion coefficient with roughness, one can show that the unfolding 
kinetics for a two-state folder deviates substantially from an Arrhenius behavior as follows. 

\ogk(f,T) = a + b/T-e 2 /T 2 (13) 

where e 2 is a constant even if the coefficients a and b change under different force and temperature 
conditions [81] . This relationship suggests that conducting force-clamp experiments over the 
range of temperatures identifies the roughness scale e. Here the condition e/AF' 1 should 
be ensured. 

To extract the roughness scale, e, using DFS, a series of DFS experiments should be performed 
as a function of T and rj so that reliable unfolding force distributions (P(f)) and corresponding 
/* value are obtained. 

k B T k B T Ax* e 2 

One way of obtaining the e from experimental data is as follows [82]. From the /* vs logr/ 
curves at two different temperatures, T\ and T 2 , one can obtain 77 (Ti) and rf(T 2 ) for which the 
/* values are identical. By equating the right-hand side of the expression in Eq|5] at T\ and T 2 
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the scale e can be estimated (EH E2] as 
2 Aa;t(Ti)fc B ri x ^{T 2 )k B T 2 
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(15) 



This equation has been used to measure e for GTPase Ran— Importin /3 complex (e > 5ksT) 
[82] and transmembrane helices (e ~ 4 — 6 fcsT) [83J. 



Pulling Speed Dependent Unfolding Pathway 

Dynamics of polymer is extremely intricate due to multiply entangled length and time scales. 
Polymeric nature of RNA and proteins immediately lends itself when a molecule interacts with 
an external force. A biopolymer adapts its configuration in response to an external stress in a 
finite amount of relaxation time. In both AFM and LOT experiments, a force is applied to one 
end of the chain with the other end being fixed. A finite amount of delay is expected for the 
tension / to propagate along the backbone of molecule and through the network of contacts that 
stabilize the native topology. To understand the effect of finite propagation time of the tension 
on unfolding dynamics, it is useful to consider a ratio between the loading rate Vj and the rate 
at which the applied force propagates along polymer chain ty (A = r^/rj). is controlled by 
experiments; r T most likely depends on the topology of a molecule. Depending the value of 
parameter A, the history of dynamics can be altered qualitatively. If A 3> 1, then the applied 
tension at one end propagates rapidly so that, even prior to the realization of the first rip, force 
along the chain is uniform. In the opposite limit, A<1, the tension is nonuniformly distributed 
along the backbone at the moment any of rupture event occurs (see the gradient of red color in 
the one at the highest loading rate, the top panel of Fig.2A). In such a situation, unraveling of 
RNA begins from a region where the value of local force exceeds the tertiary interactions. 

The unfolding simulation using SOP model of Azoarcus ribozyme provides a great insight into 
the issue of force propagation and r/-dependent unfolding pathways [67] . The intuitive argument 
given in Fig.2A is clarified by visualizing the change in the pattern of force propagation for 
Azoarcus ribozyme under three different loading conditions. Alignment of the angles between 
the bond segment vector (r i i+1 ) along force direction can tell the magnitude of force exerted at 
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each position along the backbone. The nonuniformity in the local segmental alignment is most 
evident at the highest loading rate. The dynamics of the force propagation occurs sequentially 
from one end of the chain to the other at high r/. The alignment of segment along / gets 
more homogeneous at lower Tf. These results highlight an important prediction closely related 
to polymer dynamics, that the unfolding pathways can drastically change depending on the 
loading rate, rf. By varying A (i.e., controlling rf) from A < 1 to A > 1, force experiments 
will show the dramatic effect of pulling speed dependence on the unfolding dynamics. There 
may be a dramatic change in unfolding mechanism for two different instruments using dinstinct 
different rj (LOT and AFM experiments). In addition, predictions of mechanism for forced 
unfolding based on all-atom MD simulations [M] should also be treated with caution unless due 
to topological reason, the unfolding pathways are robust to large variations in the loading rates 
regardless of A value. 



Effect of Molecular Handles on the Measurement of Hopping Transition Dynamics 

While the idea of SM experiment is to probe the dynamics of an isolated molecule, noise 
or interference from many possible sources is always a difficult problem to deal with in nano- 
scale measurement. To accurately measure the dynamics of a test molecule the interference 
between molecule and instrument should be minimal. In optical tweezers experiment, dsDNA 
or DNA-RNA hybrid handles are inserted between the microbead and test molecule so as to 
minimize the systematic error due to the microbead-molecule interaction. Unlike force-ramp, 
issues concerning non-equilibrium relaxation such as the delay of signal or force propagation 
does not lend themselves in force clamp experiments as long as the time scale of molecular 
hopping (rhop) at given force is greater than the relaxation times of handle (r^), microbead (r&), 
and other part of the instruments (typically Th op 3> T h ,r b ). Force-clamp experiment essentially 
creates an equilibrium condition. As a result, tension is uniformly distributed over the handle 
as well as the test molecule while the molecule hops. 

However, critical issue as to measurements still remains even at a perfect equilibrium due 
to handle fluctuations, (i) Since the dynamics of test molecule is measured by monitoring the 
position of microbead, to gain an identical signal with the test molecule, the fluctuation of the 
handle should be minimal. Therefore, one may conclude that short and stiff handles (smaller 
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L/l p ) are ideal for precisely sampling the conformational statistics j3l] (see Fig.3B and caption), 
(ii) Simulations of hopping dynamics with handles differing length and flexibility shows that the 
hopping kinetics is least compromised from the true handle-free kinetics when short and flexible 
handles are used (85]. Physically, when a test molecule is sandwiched between handles, the 
diffusive motion is dynamically pinned. The stiffer handles, the slower the hopping transition. 
Therefore, short and flexible handles are suitable to preserve the true dynamics of the test 
molecule. 

The above two conditions for (i) precise measurement of thermodynamics and (ii) accurate 
measurement of true kinetics apparently contradict each other. However, this dilemma can be 
avoided by simply measuring the folding landscape accurately using short and stiff handle as 
long as the effective diffusion coefficient associated with the reaction coordinate is known. At 
least, at / ~ f m id, the hopping time trace can easily fulfill the ergodic condition, providing a good 
equilibrium free energy profile. Furthermore, it can be shown that the molecular extension z (or 
the end-to-end distance R) used to represent the free energy profile is indeed a good reaction 
coordinate on which the Bell's prescription to obtain rate agrees well with the folding rate one can 
obtain using a simulation; for a given transition rate from NBA to UBA kpifmid), the variation 
of force from / = f mid modifies the rate as k F (f mid ± 5f) = k F (f mid ) exp (±5f ■ Ax t /k B T). 
Once an accurate free energy (F(z m ) ~ F°(z m ) ps F(z sys ) holds for L/l p <C 1. see Fig.3B and 
caption) is obtained, one can determine the hopping kinetics by directly calculating the mean 
first passage time on F°(z m ). The only unknown, the effective diffusion coefficient on F°(z m ), 
can be estimated by mapping the hopping dynamics of test molecule to the dynamics of a simple 
analytical model such as generalized Rouse model (GRM) [55"] 157]. 

Folding Dynamics upon Force-Quench 

Because of the multi-dimensional nature, the dynamics of biomolecules is sensitive to the 
condition to which the molecule is imposed. In a rugged energy landscape, the folding rate is 
not unique. Rather, the folding kinetics can vary greatly depending on initial conditions [13J. 
The difference in the initial condition may be due to either urea concentration, temperature, or 
force. For RNA it is well known that the initial counter ion condition can alter the folding route 
in a drastic fashion [13J. A biopolymer of interest adapts its structure to a given condition and 
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creates a "condition specific denatured state ensemble (DSE)." The routes to the NBA from the 
DSE are determined by the shape and ruggedness of energy landscape due to the factors such 
as side-chain interactions, topological frustrations and so forth. 

The force quench refolding dynamics of poly-ubiquitin (poly-Ub) by Fernandez and Li [11] 
is the first experiment that focused on the folding dynamics of proteins from a fully extended 
ensemble, in which the dynamics of molecular extension of poly-Ub construct was traced from an 
initial stretched ensemble prepared at high stretch-force fs = 122 pN to a quenched ensemble 
at low quench- force fg = 15 pN (see Fig.4A). The folding trajectories monitored using the 
molecular extension (R) was at least an order of magnitude slower than the one from bulk 
measurement and was characterized with continuous transitions divided into at least three (four 
in [UJ) stages: initial reduction of R, long flat plateau, and a cooperative collapse transition at 
the final stage. Two main questions were immediately raised from the results of experiments, 
(i) Why is the folding (collapse) process so slow compared to the one at bulk measurement? (ii) 
What is the nature of the plateau and cooperative transition at the final stage? Concerning the 
point (i), it was suspected that an aggregate was formed between the Ub monomer since the 
effective concentration is more than the critical concentration for the aggregate formation [88J . 
But, this possibility was ruled out since no signature of disrupting the aggregate contacts was 
found when the folded poly-Ub construct was re-stretched; the number of step indicating the 
unfolding of individual Ub domain was consistent with the number of Ub monomers [H] . 

The immediate interpretation toward the anomalous behavior of refolding trajectories upon 
force quench is found from the vastly different initial structural ensemble [12] H9l ETJ [89] [90] . 
The initial structural ensemble under high tension is fully stretched (SSE, stretched state en- 
semble) while the nature of ensemble for bulk measurements are thermally denatured (thermally 
denatured ensemble, TDE). In terms of structural characteristics, TDE and SSE have drastic 
difference. The entropy of SSE is smaller than TDE (Ssse < Stde) P21 E7J [90]. Therefore, it 
is not unusual that the folding kinetics upon force quench is vastly different from the one at the 
bulk measurement. It is conceivable that the energy landscape explored starting from these two 
distinct ensembles vastly differ. 

The next more elaborate question is then why the folding rate for force-quench is slower. 
Under force quench condition, two driving forces compete each other. One is the free energy 
gradient bias toward the NBA, the other is the quench force (/q < f m id) that resists the collapse 
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process. A free energy barrier is formed under these two competing forces. Consequently, prior 
to making a transition to the NBA, the molecule is trapped in a finite sized free energy barrier, 
forming a metastable intermediate. Due to fq, the formation of contact responsible for the 
collapse process is suppressed. This tendency increases with an increasing fq, which again 
increases the refolding time (t f ). If the time spent for force quench (tq) is too fast, i.e., 
tq <C Tp, producing a nonequilibrated system in which the molecule is trapped in the force- 
induced metastable intermediate (FIMI), then the plateau in terms of molecular extension x 
will be observed in the measurement. Compared to TDE where the contacts responsible for 
the collapse are in proximity, the formation of folding (or collapse) nuclei is much more time- 
consuming for a system trapped in FIMI state at higher fq. Thus, the folding route from SSE 
to NBA differs significantly from the one from TDE to NBA [90J. 

The FIMI is generic to the force-quench refolding dynamics of any biopolymer, which can 
be easily tested by applying a tension to a semiflexible polymer that forms a toroid or racquets 
structures upon collapse. The energy hamiltonian is given by 



with the parameters, eu = 1.5ksT, k s = 2000/csT, iV = 200, a = 0.6 and kb = 80 k^T . By 
abruptly changing / = f s = 83 pN to f — fq — 4 pN, long plateaux with varying durations 
are observed in the semiflexible polymer described by the above energy hamiltonian (Fig.4B). 
An optical tweezers experiment showing the difficulty of forming DNA toroid under tension 
unambiguously supports the argument that FIMI is generic to collapse dynamics of any polymer 
under tension [91] (Fig.4C). 

Concluding Remarks 

The effort to decipher the energy landscape of a biopolymer by monitoring its mechanical 
response to an external stress has greatly enriched single molecule experiments, theory and 
molecular simulations associated with force mechanics of biomolecules and cells. We have also 
witnessed many instances of successful application of theories and molecular simulations for the 
experimental data analysis [EH EH E21 [92]. However, as experiments are conducted on more 
complicated and larger systems among cellular constituents, non-trivial patterns of molecular 
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response are more observed [621 I93H95] . Theoretical method to analyze the mechanical response 
from multiple barriers and studies on the structural origin of more non-trivial mechanical 
responses such as catch-slip bond are still on their infancy. Further studies need to be done 
beyond force induced dynamics of a molecular system associated with a single free energy 
barrier crossing dynamics on a one- dimensional energy profile [HI W2\ l9"6"] . Interference between 
a molecule of interest and instrument itself pose the problem of deconvolution [Ml [55]. Of great 
interest from the perspective of both theory and simulation will be revealing how hydrodynamic 
interaction between the subdomains of a biopolymer with complex architecture play a role 
during its forced-unfolding process. More exciting experiments using mechanical force and 
breakthroughs in theories and molecular simulation methods are anticipated in the next decade 
to further reveal the beauty of living systems at the microscopic level. 

I am grateful to Dave Thirumalai for intellectually pleasurable discussions and collaboration 
on theory for force spectroscopy and molecular simulations for the last many years. This work 
was supported in part by grants from National Research Foundation of Korea (RO 1-2008-000- 
10920-0). 
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Figure Captions 

Fig.l Force-ramp unfolding simulation of an RNA hairpin using SOP model. (A) FEC 
(/(-R)) at rf = 45pN/s (k = 0.07pN/nm, v = 0.64/im/s). f(R) can be decomposed using time 
t into R(t) and /(£). The experimental setup of optical tweezer is mimicked by attaching a 
harmonic spring with strength k. The force is recorded by measuring the extension of this 
spring. f(t) shows how the force is ramped with time when the system is pulled with constant 
pulling speed v. The two-state hopping transition in R and / begins when / reaches 12 
pN and ends at / ~ 17 pN. (B) Unfolding force distribution from 100 trajectories for each 
77. rj = r° f , r° f /3, 20r?/3, 10r° f when r° f = 4.5 x 10 3 pN/s. (C) [logr/, /*] plot obtained from 
(B), which clearly shows a positive curvature. Linear regression at different rj values provide 
different Ax$ values. Figure adapted from pjg] . 

Fig. 2 Force propagation. (A) A diagram depicting the pulling speed dependent unfolding 
pathway. For a biopolymer consisting o two hairpins 1 and 2 whose barriers associated with 
forced unfolding are given AFl and AFj with AF} < AF|, the unfolding will occur through one 
of the two reaction paths depending on the rf. For A > 1, tension will uniformly distribute along 
the chain, so that fi w f 2 . Since AP* — f\Ax$ < APj — f 2 Ax$, the hairpin 1 will unfold before 
hairpin 2 (pathway I). In contrast, for A <C 1, /i < /2 leads to AP* — f%Ax* > AF 2 — f 2 Ax*, 
thus the hairpin 2 will unfold before hairpin 1 (pathway II). (B) rj-dependent unfolding 
pathway simulated with Azoarcus ribozyme using SOP model. FEC (middle panel) shows that 
unfolding occurs via N -> [P5] -> [P6] ->■ [P2] -> [PA] -> [P3] -> [PI] at r/ = 1.2 x 10 6 
pN/s (red). -»■ [P1,P5,P6] -»■ [P2] [P4] [P3] at r/ = 3.6 x 10 5 pN/s (green), and 
N -> [P1,P2,P5,P6] ->■ [P3,P4] at r/ = 1.8 x 10 4 pN/s (blue). From high to low r f the 
unfolding pathways were changed completely. Time evolutions of cos#i (z = 1, 2, . . . , N — 1) at 
three r/ visualize how the tension is being propagated. Figure adapted from [67] 

Fig. 3 Handle effect on the measurement and RNA hopping dynamics. (A) Molecular 
simulation of RNA hopping dynamics with two handles of L = 25 nm and l p = 70 nm attached 
to the 5' and 3' ends. The illustration was created by using the simulated structures collected 
every 0.5 ms. An example of the time trace of each component of the system, at / = 15.4 pN. 



27 



z m and z sys measure the extension dynamics of the RNA hairpin and of the handle- RNA-handle 
system, respectively The time-averaged value z r (t) = (1/r) J^dTz sys {T) for the time trace 
of z sys is shown as the bold line. The histograms of the extension are shown on top of each 
column. (B) The free energy profiles, F eq (z sys ) from the two ends of handles attached (dashed 
line in blue), F eq (z m ) from 5' and 3' ends of RNA (solid line in red), and F° q (z m ) from 5' and 
3' ends of RNA with no handles attached. When L/l p <gC 1, F eq (z S y S ) ~ F eq (z m ) « F°(z m ). 
(C) With increasing handle length the folding rate (k'p(L) deviates from a true folding rate 
(kp(0)). This effect is larger for handles with a greater l p . Figure adapted from [86] . 

Fig. 4 Refolding (Collapse) dynamics of biopolymers upon force quench. (A) Refolding 
trajectory of poly-Ub generated by atomic force microscopy upon force quench (Figure adapted 
from pT|). (B) The collapse of a semiflexible chain in a poor solvent under force quench 
condition. Toroidal (magenta), racquet structures (green) are formed after varying durations of 
plateaux associated with FIMIs (Figure adapted from [12J). (C) After an addition of multivalent 
counterions, it takes as long as ~ 2400 sec for a A-DNA to initiate collapse dynamics under a 
small quench force /q = 1.2 pN (Figure adapted from pTj). 
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